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Abstract 



Contrary to what is often stated, a fundamental spacetime discreteness need 
, not contradict Lorentz invariance. A causal set's discreteness is in fact locally 

lO I Lorentz invariant, and we recall the reasons why. For illustration, we introduce 

■ a phenomenological model of massive particles propagating in a Minkowski 

spacetime which arises from an underlying causal set. The particles undergo 
a Lorentz invariant diffusion in phase space, and we speculate on whether this 
' could have any bearing on the origin of high energy cosmic rays. 

7' 

u 

In discrete approaches to quantum gravity, the fundamental description of space- 
time is not taken to be a manifold, but some discrete structure to which the manifold 
^ I is only an approximation. The scale of this discreteness is usually assumed to be 

I Planckian. It is often asserted any such theory must violate local Lorentz invariance 

(LLI) and a new area of research - LLI violating phenomenological effects of quantum 
gravity - has grown up around this idea. The purpose of this letter is to emphasize 
that causal set theory P respects LLI and to open a new phenomenological window 
on this approach to quantum gravity. 

What does it mean to say that a discrete theory respects Lorentz invariance? It 
is difficult to give a precise answer, but intuitively the import is clear. Whenever a 
continuum is a good approximation to the underlying structure (and assuming specif- 
ically that the approximating continuum is a Lorentzian manifold M), the underlying 
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discreteness must not, in and of itself, suffice to distinguish a local Lorentz frame at 
any point of M. In consequence, no phenomenological theory in M derived from such 
a scheme can involve a local (or global) Lorentz frame either.^ 

Of course the above presupposes an answer to the question: "How is the approx- 
imating continuum related to the discrete entity that underlies it?". Whether or 
not a particular discrete theory respects LLI cannot be settled until this question is 
answered in the context of that theory. Luckily, in causal set theory, there is a clear 
proposal for an answer, and we will show that LLI is indeed respected. 

A causal set (causet for short) is a locally finite, partially ordered set (for reviews 
and motivation for causal set theory see |2] |S1)- This is a set, C, endowed with 
a binary relation -< such that elements of the set satisfy the conditions (i) (x -< 
y)a.iad{y -< z) =^ {x -< z) (transitivity), (ii) x -/{ x (acyclicity), and (iii) all 
intervals {x : y -< x -< z} are finite. The relation -< gives rise, in the continuum limit, 
to the causal order on spacetime points, and the number of elements in a subcauset 
yields the volume of the corresponding region of the continuum in Planck units. 

In the continuum context, the causal order and volume information suffice to 
specify a (causally reasonable) Lorentzian manifold jH E] . It is therefore reasonable 
to regard a Lorentzian manifold as an approximation to a causet if that causet is 
a discrete "sampling" of the continuum causal order with uniform density.^ More 
specifically, we may say that a Lorentzian manifold M approximates a causet C 
(M !^ C) if C could have arisen, with relatively high probability, via a random 
process of "sprinkling into M" , at Planck density,^ with the causet relations induced 
by the spacetime causal structure.^ 

A "sprinkling" is more properly described as a Poisson process. To see what this 
means, imagine dividing M, using any local coordinate systems, into small boxes of 
volume V , and then placing a "sprinkled point" independently into each box with 
probability V/Vfundi where Vfund is the fundamental volume (of order the Planck 
volume). The Poisson process is the limit of this procedure as V tends to zero. 
Because spacetime volume is an invariant, the limiting process is independent of the 
coordinate systems used to define the boxes. It follows that one cannot tell which 
frame was used to produce the sprinkling: the approximation is "equally good in all 

^Naturally, there can be no question of a literal action of the entire Lorentz group on an individual 
discrete structure. Rather such a structure can only be Lorentz invariant in the same sense that a 
fluid is translation invariant. This should not detract from the fact that a fluid is indeed translation 
invariant in an important sense, whereas a crystalline solid is not. 

^More generally, one would only require that some coarse-graining of the causet approximate M 
in this sense; but we ignore this distinction here. For a somewhat different approach to defining a 
relation of closeness between a manifold and a causet see |S] . 

•^By "Planck density" , we really mean "density unity in fundamental units" . One expects funda- 
mental units to be equal in order of magnitude to Planck units. 

"^Taking the order relation of the causet to be induced strictly from that of the spacetime is only 
the simplest possibility. Other rules could be considered, but they would not affect anything in this 
paper. 
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frames." 




Figure 1: A regular lattice of spacetime points in two different Lorentz frames. Normal 
conventions for spacetime axes are used. While in (a) the lattice appears to have a regular 
density of elements, in the boosted frame shown in (b) the density of points is revealed not 
to be uniform. 

Why is the randomness crucial? Let us take the example of 1+1 dimensional 
Minkowski space. One obvious way to try to discretize it is to choose a frame and 
use a "diamond lattice" with respect to that frame, i.e. the points with coordinates 
{t,x) = (e(r + s),e(r — s)), where r and s are integers and e is some fixed length, 
as in figure (^). In this frame the lattice appears to be a good approximation; all 
"nicely shaped" large regions have a similar density of elements. In a frame boosted 
at velocity V in the positive direction, however, the elements are at (7e[(l — V)r + 
(1 + V^)s],7e[(l - V)r - (1 + V)s]). Figure jib) shows this lattice with 7 = 1.25. 
Now it becomes clear that, if the boost is large enough, there will be "nice" big 
regions containing no elements at all, and others containing far too many elements. 
The approximation only looks good in the original frame, and so it breaks Lorentz 
invariance by preferring this frame. In light of this example, it seems likely that the 
same problem would affect any non-random discretization of a spacetime. Thus, for 
example, a Regge-type triangulation whose simplices look "fat" in one frame will look 
"long and skinny" in a relatively highly boosted frame. 

We want to emphasise that not only is the process of sprinkling Lorentz invariant 
but so also are almost all of the individual causets that are generated. An objection 
that often comes up in this connection concerns the necessary occurrence of voids 
in any given Poisson sprinkling. While it is true that voids must occur, this does 
not cause a problem for Lorentz invariance (or any other problem that we know of). 
However, one may still feel uneasy about the voids, and some people seem to believe 
that they necessarily would break Lorentz invariance in some manner. To put such 
qualms in perspective, let us estimate the probability that there is at least one void 
of nuclear dimensions in the history of the observable universe since the Big Bang. 
More precisely, we will bound the probability that a sprinkling would leave empty 
any interval whose height is of the order of one Fermi. (An interval in spacetime is a 
"double light cone" or "Alexandrov neighborhood". We do not require the "axis" of 
the interval to be aligned with the cosmic rest frame. Hence our bound will apply to 
the probability of "finding a void in any frame".) 

All the numbers in what follows are "of the order of" . Consider as a model of the 
universe a portion P of Minkowski spacetime, the size of the observable universe and 
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defined by < t < T, < < T, z = 1, 2, 3 in some frame. If T is 13 billion years, 
the spacetime volume is 10^'*'' in fundamental units. An interval of nuclear size has 
spacetime volume 10'^''. This means that the probability that any particular nuclear 
sized interval will be a void is e~^° . But we want the probability q that at least 
one interval (any one) in P will be void. This can't be calculated easily, because 
the intervals overlap and the probabilities for them to be void are not independent. 
However, we can put an upper bound on q without much difficulty. 

Let us fill P with coordinate balls (with respect to the "defining frame" of P) 
of small enough radius that any "upright" interval of nuclear size is guaranteed to 
contain at least one complete ball. (Upright means that the top and bottom points 
have the same spatial coordinates.) A radius of one hundredth nuclear size will do, 
and we will have 10^^^ of these balls packed into P. The probability that at least one 
of them is void is 10^^^e~^° . These same balls will also suffice for intervals of nuclear 
volume that are slightly boosted or "tilted" from the upright. They will certainly do 
for all 7 factors less than or equal to 7 = 5/4, that is for a region of the Lorentz group 
with volume of order 1. For each such cell of the Lorentz group we choose a set of 
coordinate balls in spacetime (in the corresponding frame) . The relevant region of the 
Lorentz group is bounded by the maximum relevant boost in P, which corresponds 
to a 7 factor of 10''^. (Any larger boost would produce a ball that could not fit into 
P. In fact the maximum 7 is actually smaller, since the boosted nucleus would meet 
the boundary of P before the smaller ball would.) The number of cells needed to 
cover this region of the Lorentz group is 10^^. The probability of getting any nuclear 
sized void is less than the probability that any one of the coordinate balls from any of 
the boosted sets will be void. This in turn is less than 10^^ x 10^^^xe~^° , a number 
so tiny that the two prefactors have no impact whatsoever on its value. 

Given that causal sets respect Lorentz invariance, what conclusions can be drawn? 
Most obviously, we predict that no violation of LLI will be observed at the phe- 
nomenological level, so that, if any of the experiments currently planned or under 
way did find such a violation, the causal set hypothesis would be disfavored. But this 
is only a negative prediction. Are there also positive signatures? Since the causal set 
hypothesis makes such a definite statement about the underlying structure for space- 
time and how it is related to the continuum we actually experience, it is not difficult 
to devise concrete models predicting potentially observable effects of the underlying 
discreteness. We give one such model below which we call "swerving" , after Lucretius 

m 

"The atoms must a little swerve at times - but only the least, lest we should seem 
to feign motions oblique, and fact refute us there." 

In the continuum, massive particles travel on timelike geodesies. However, an 
underlying discreteness might induce small fluctuations in the particle's worldline, 
and the causal set picture naturally suggests models in which this effect would be 
Lorentz invariant. Though it might be too small to observe on everyday scales, such 
an effect might be detectable in sensitive laboratory experiments, or by astronomical 
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observations if the particle were travelling over cosmic distances. Here is a model of 
this type. 

Consider a hypothetical point-particle of mass m moving through a causet C 
derived by sprinkling Minkowski space M^. We will take its trajectory to be a chain 
of elements of C (i.e. a totally ordered subset of C). No such chain can correspond 
perfectly to a straight line in M^. So, given an "initial segment" of the trajectory 
up to some element e„ of C, how could its future continuation be determined? If we 
assume that the trajectory's past determines its future, but that only a certain proper 
time Tf into the past of the trajectory is relevant, then we are led to the following 
Lorentz invariant rule as a particularly simple discrete analog of geodesic motion. 




Figure 2: A portion of 1 + 1 Minkowski space wtiere tlie dots represent elements of a causet 
sprinkled into it. The trajectory of the particle has reached with momentum pn (the 
frame having been chosen so that the three-momentum is zero). The dotted line is the 
hyperbola of points a proper time tj to the future of e„. The element within proper time 
Tf of e„ that best preserves the momentum is e^+i. The ratio of tJ to the density of 
sprinkling has been exaggerated here to emphasize the momentum change in one step. In 
a more realistic model Tf would be larger. 

The future trajectory is constructed inductively, as illustrated in fig. |21 Starting 
from an element e„ and a momentum p^, the next element in the trajectory, e„+i, 
must be chosen. For convenience we have drawn the diagram in the rest frame at e^. 
The new momentum Pn+i is defined to be proportional to the vector between and 
Cn+i- The selection of e„+i is made such that e„+i is in the causal future of e„ and is 
within proper time Tf of and so that \pn+i ~ Pn\ is minimized. In the figure, this 
means that e„+i is within the future light cone of e„, below the dotted hyperbola, and 
such that the vector from e„ to e„+i is as close to the vertical as possible. (That there 
exists such an element in that region is guaranteed by the infinite spacetime volume of 
the region, the Poisson distribution and the local finiteness condition.) This realizes 
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the ideas that the trajectory should be as close to a straight line as possible and that 
the dynamics should be approximately Markovian if the "forgetting time" tj is small. 
(There is ample room for it to be small and yet much bigger than the discreteness 
scale if the latter is Planckian.) The process is then repeated starting with e„+i and 

Pn+l- 

Our model implies random fluctuations in the momentum of the particle. For 
any numbers 5i > ^2 > at stage n there is a finite volume within proper time Tf 
to the future of e„ such that, if element e„+i were in this volume, the momentum 
change \pn+i — Pn\ would lie between 6i and 62. The probability of this happening is 
the probability that this volume is not empty of sprinkled points, while the volume 
leading to a smaller momentum change is empty. Both these probabilities are given 
precisely by the Poisson distribution: the probability that a volume U is empty is 
e~^^, where g is the density of the sprinkling. 

This simple model can be criticised on many grounds: it treats the particle as if 
it were of zero size, it is deterministic rather than quantum, etc.. Most seriously, it 
could not possibly be fundamental, since the law of motion of the trajectory is not 
formulated in terms of the causet, but refers also to the approximating Minkowski 
spacetime.^ We present it in the spirit of |l8j: form a concrete model with testable 
consequences based on important aspects of the fundamental formalism and compare 
to observation. 

Since the hypothesized Lucretius effect is supposed to occur on very small scales, 
it should be possible to approximate it over macroscopic distances by a diffusion 
equation (hydrodynamic limit). This is analogous to how the ordinary diffusion 
equation describes the long time behaviour of a random walk. In our case however, 
the diffusion is not in physical space R^, but on the phase space x M^, where 
is the mass shell (Lobachevskii space), is Minkowski spacetime, and the diffusion 
takes place in proper time t. The diffusion in spacetime is secondary and is driven 
by that in momentum space, in close analogy with the Ornstein-Uhlenbeck process. 

Consider the scalar (not scalar density) probability distribution p = p{p'^ ,x^]t) 
on X M^. It is a function of momentum, p", spacetime position, x^, and proper time 
T. We write the full four-momentum as an argument with the understanding that 
on the mass shell there are only three independent components. With the condition 
that the process be Lorentz invariant, the following equation for p can be derived 
by following the prescription set out in for stochastic evolution on a manifold of 
states: 

where Vp is the Laplacian on tf, m is the mass of the particle, and is a constant 
(that will depend on the parameters of the discrete process, such as the forgetting 
time Tf). 

^This defect can be overcome fairly easily, however. 
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This equation defines a diffusion process in which the particle's proper time r 
serves as time. It is the unique Markovian, Poincare invariant, relativistically causal, 
diffusion law in vacuum that preserves p > 0.^ For fixed mass, it has a single free 
parameter, the diffusion constant k. (In principle, the coefficient of the P^-^p term 
could be different, but that would mean that would not be the physical mo- 
mentum.) Given these uniqueness properties, the equation should be insensitive to 
variations in the microscopic model that underlies it, so long as the latter is Lorentz 
invariant, causal and (approximately) Markovian. So although the above discrete 
swerve model might be wrong in detail, the macroscopic phenomenology of (Q) tran- 
scends it. 

If p is a function of momentum alone and is initially (at r = 0) a delta function 
in momentum, then the (un-normalized) solution - adapted from the solution of the 
diffusion equation on pT, - is 

p(j9) = e-^'/4^^(A;r)-ie-^"— ^ (2) 

where k = k/m^c^, R = sinh~^(p/mc), and p = |p| is the norm of the three- 
momentum in the frame defined by the point in at which the diffusion begins. 
a = mcR is then the geodesic distance from that point. 

Equation ((T)) and its fundamental solution Q, expressed as they are in terms 
of proper time, are not well suited to comparison with experiment/observation, even 
though they exhibit the underlying Lorentz invariance very clearly. Instead, one needs 
a description of the same process with respect to cosmic or laboratory time. To this 
end, fix a preferred set of spacelike hypersurfaces t = = constant, and assume that 
the distribution on some initial hypersurface of the set is uniform in space {i.e. p 
does not depend on x\ where i are spatial coordinates). The swerves, being spatially 
homogeneous and isotropic, will preserve this uniformity, and so if we assume that 
any additional frictional effects are also homogeneous and isotropic, the distribution 
will remain uniform. We require, under these conditions, an equation governing how 
the probability distribution evolves in cosmic time t. 

Such an equation can be deduced from ([T]). (Details of the derivation and of the 
inhomogeneous case will appear elsewhere.) The result is: 

^ = kVl( , ^ 1 - V,(^» (3) 

Here, the scalar function p = p(p^; t) on tf, gives the momentum distribution, Va is 
the covariant derivative on H^, and p is the norm of the particle's three-momentum, 
p in the cosmic frame. (To make this equation plausible, notice that the factor 

^The Ornstein-Uhlenbeck process for a particle diffusing in interaction with a relativistic fluid 
has been considered in jlO| . In that case the rest frame of the fluid provides a preferred frame for 
the stochastic noise term driving the diffusion. The present process has no preferred frame and is 
fundamentaUy Lorentz invariant. 
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\/l + /w? is the boost factor 7 = dt/dr.) The term involving the vector w°' is 
a friction term added to represent the effect on the particle's momentum of, for 
example, the Hubble expansion and interactions with the CMBR (cosmic microwave 
background radiation). The specific form of w, which will in general be a function of 
the momentum, will depend on the type of friction involved/ For the cases mentioned, 
w"" will have only a radial component (in the p = |p| direction). 

In the time- dependent case, Q will probably have to be solved numerically. How- 
ever, we can hope to analyze its equilibrium solutions analytically. For example, at 
large p, energy E p, and if ~ —bE"' with n > 1, then the equilibrium solu- 
tion will behave as Ee~^^"^^. If the dominant friction over any high-energy range is 
constant, i.e. if w'^ ~ dE/dt ~ constant, then the equilibrium distribution will be a 
power law in that range. 

How might one observe diffusion of the above sort? Cosmological and astro- 
physical observations are the obvious places to look for consequences of a universal 
acceleration mechanism. But first, laboratory physics can put an upper bound on the 
diffusion constant k. Suppose the particles in question are protons. If k were large 
enough, hydrogen gas would spontaneously heat up in a short time, and this has not 
been observed. In the laboratory regime, hydrogen is non-relativistic, so equation Q 
can be approximated as 

I - >=^'^ w 

where is now the standard Laplacian on R^. This is the standard diffusion equation 
and has the well known solution: 

p = A(t)exp(-^) (5) 

where A{t) is a normalization factor. Usefully, this is also the form of the Maxwell 
distribution for a classical gas in thermal equilibrium: 

PMaxweU = Aexp(-^-£-^) (6) 

where m is the molecular mass, T is the temperature, and ks is Boltzmann's constant. 
If a gas starts in a thermal state, it will therefore remain in a thermal state even if 
swerves are included. Moreover, the above two equations imply that the temperature 
will scale linearly with time, specifically: 

dT 2k 

-r = —r- (7) 

dt mkB 

Assuming, for the sake of argument, that a heating rate of a millionth of a degree per 
second would already have been detected in the laboratory, we obtain the approximate 
bound 

k < IQ-^'^ kg'^mh-^ (8) 

'^In the case of violent momentum transfers, this would have to be generalized to a Boltzmann 
type collision term. 
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The maximum average energy gain due to swerves consistent with this rate of tem- 
perature gain can be obtained from the formula {E) = 3A;bT/2: 

{AE)/At < 4.3 X 10-^^eVs-^ (9) 

Now, let us turn to some possible astrophysical effects of swerves. One outstanding 
astronomical puzzle is the origin of high energy cosmic rays (see ^2] for a recent 
review). Attention is often focused on the so-called "trans-GKZ" events, apparent 
detections of cosmic rays with energies above 5 x 10^^ eV. Such primaries, if they 
are protons, cannot have come from farther than about 20 Mpc (because they would 
have decayed due to photo-pion production with the CMB photons), but they have 
no obvious source in that distance range. But even for cosmic rays between 10^^ 
eV and 10^^ eV, there are only suggestions and no universally accepted acceleration 
mechanism for producing the observed energy distribution. The data (see e.g. fig 1 of 
|12j ) seem to cry out for a universal cosmic acceleration mechanism that would inject 
protons, say, into the galaxy with a power law distribution of where 2 < a < 3, 
so that the observed variations in the power law and deviations from isotropy would 
be due to the dynamics of the protons in the galaxy. Could swerves provide such a 
cosmic mechanism?^ 

Swerves induce a "statistical acceleration" analogous to Fermi acceleration, and 
it is possible a priori that enough intergalactic hydrogen could be accelerated up to 
very high energies to explain the data. Unfortunately this degree of acceleration is 
inconsistent with the bound on k already discussed. A "Lucretian" explanation of 
the cosmic ray data, assuming the primaries are protons, would require some protons 
to accelerate to ~ 10^° eV from far lower energies on a timescale of the age of the 
universe. To produce a power law distribution in energy, a significant proportion of 
those protons reaching, say, 10^^ eV would have to go on to double their energy and 
more. The rest frame of a proton with an energy of 10^^ eV has a 7 factor of 10^ 
relative to the cosmic frame, so 10 billion years of cosmic time is only 10 years of 
proper time for such a proton. Doubling its energy in the cosmic frame would mean 
gaining around 250 MeV in its own frame. But from the inequality in this frame 
the average energy gain in 10 years could be at most 1.4 x 10~^ eV. (At these energies, 
we can trust the non-relativistic approximation.) Since the distribution of momentum 
is Gaussian, with such a low average energy gain, the probability of gaining 250 MeV 
is exponentially small. In other words, a proton has practically no chance of making 
it from 10^^ eV to 2 x 10^*^ eV in the age of the universe as a result of swerves. This 
calculation assumes that k is roughly the same for an intergalactic hydrogen atom as 
it is for a proton, as it is for a H2 molecule in a box of gas, but the argument is still 
valid even for a k many orders of magnitude larger than has been assumed. 

Proton swerves cannot explain trans-GKZ cosmic rays either (if indeed their ap- 
parent observation turns out to be correct). Swerves would accelerate some protons 

^The idea that the rays might be the result of spontaneous acceleration, as a result of non- 
standard QFT, has been discussed in 
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back up beyond 10^° eV after they entered the "GKZ sphere" , ~ 20 Mpc from us; but 
this effect would not be significant. The argument for this is similar to that above. 
A proton with an energy of 10^^ eV would reach us from the GKZ sphere in about 2 
hours of proper time. This would not leave enough time for a non-negligible fraction 
of protons to, say, double their energy. 

So the most direct application of the swerve idea to protons cannot explain the 
origin of high energy cosmic rays. However, more complicated scenarios can be 
considered. For example, in [T3j the authors postulate homogeneously distributed 
sources producing (by some unknown mechanism) neutrinos with energies above 10^^ 
eV which collide with a cosmological background of neutrinos (hot dark matter) to 
produce - amongst other things - protons and gamma rays that could be cosmic 
ray primaries. Perhaps swerves could provide the required acceleration in this case. 
Indeed neutrinos are more likely candidates than protons to be affected by the un- 
derlying discreteness as they are more point-like than protons, according to present 
beliefs. Moreover, we have few if any laboratory bounds on k to contend with in the 
case of neutrinos. 

More sophisticated models could also be developed. Our simple proton model 
assumes that the "diffusion constant" k does not depend on local factors like average 
particle density, temperature etc. In a more realistic swerve model, perhaps a high 
particle density would lower the rate of diffusion, in which case the constraints from 
laboratory physics could be loosened. A second sort of generalization of the diffusion 
(or Fokker-Planck) equation Q would relax the assumption of locality in momentum 
space. The "friction term" (jS)) captures the effect of many small momentum transfers 
due to (for example) CMBR scattering. The effect of large kicks would have to be 
described by a Boltzmann equation. A third improvement to our model would be to 
treat the particles quantum mechanically rather than classically, allowing one to take 
into account the finite size of the "wave packet". Such a change might be important, 
because it is conceivable that matter-induced decoherence would influence the value 
of k, making it different on earth than in interstellar space. Unfortunately, however, 
the type of over-arching framework that is available for classical diffusion seems to 
be lacking in the quantum case, so it is less obvious how to proceed. Finally, in a full 
theory of causal set quantum gravity, regions of continuum spacetime might be best 
described as a quantum superposition of many causets, and a better phenomenological 
model might have to reflect this aspect as well. 

Let us return for a moment to possible observational evidence for swerves. Since 
the path of a particle would no longer be an exact geodesic, a certain amount of 
fuzzing of distant sources of particles would occur. Perhaps this could be revealed by 
highly directional detectors of some sort. 

So far we have limited our discussion to the case of massive particles. A Lorentz 
invariant diffusion equation for massless particles can also be written down, although 
in this case we lack a concrete model of propagation on the underlying causet that 
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could serve as motivation. ^ Rather than diffusion in H^, we have in this case diffusion 
on a "hght" cone, since the 4-momentum of a massless particle is a null vector. We will 
describe this further in a future work. Lorentz invariant diffusion on the cone cannot 
alter the direction of the momentum, but it will cause its magnitude to fluctuate, so 
that a distribution peaked at a certain energy would spread with time. Accordingly, 
one might seek evidence for this kind of diffusion in the blurring of sharply peaked 
spectral features of distant sources (such as emission and absorption lines). This is 
in contrast to what has been proposed in Lorentz-violating models, where the speed 
of a photon is presumed to vary with its energy (see e.g. [IS]). 

Finally, just as tests of Lorentz invariance push special relativity to its limits, 
tests of unitarity would allow one to push the state-vector formalism of quantum 
mechanics to its limits. It seems only reasonable therefore that simple models of the 
possible effects of non-unitarity be formulated. A non-relativistic example is given in 
|16j . Perhaps a relativistic version of this could produce a similar effect to swerves. 

To summarize, there is no reason that the assumption of an underlying spacetime 
discreteness must give rise to violations of local Lorentz invariance, because the causal 
set hypothesis does not. To illustrate the point that discreteness can nevertheless 
have observable effects, we have exhibited a Lorentz invariant momentum diffusion 
motivated by causal sets. If it is indeed the case that certain proposals for quantum 
gravity, such as loop quantum gravity and spin foams do predict violations of LLI, 
then we are in the happy situation of having a way to distinguish between different 
proposals experimentally. Be that as it may, the specificity of the models treated in 
this paper indicates that the causal set approach holds great potential for providing 
phenomenological theories of matter propagating in a discrete background. In an era 
of ever increasing sensitivity and power in cosmological observations, this potential 
to predict and detect the effects of a fundamental spacetime discreteness should be 
exploited. 

Note added: The Poincare-invariant diffusion process described above is con- 
structed rigorously in [T7] and [TH] . 

We are grateful to Pasquale Blasi, David Craig, Stuart Dowker, Seth Major, Laura 
Mersini, Angela Olinto and David Rideout for discussing these ideas with us. 

This research was partly supported by NSF grant PHY-0098488 at SU, by Air 
Force grant AFOSR at UCSD, by the Office of Research and Computing of Syracuse 
University and by the Department of Physics at Queen Mary, University of London. 

^We also lack an idea of how to describe such a diffusion in wave language as opposed to particle 
language. What would it mean in the case of Maxwell's equations for example? 
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